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1 Introduction 


The foundations for a Hamiltonian formulation of constrained systems were orig¬ 
inally laid down by Dirac |Q] whose work still remains fundamental to our under¬ 
standing of the subject. For second class systems, however, this approach poses 
problems which are related to ambiguities that may arise in the transition of (clas¬ 
sical) Dirac brackets to (quantum) commutators. During recent years a viable and 
powerful alternative for discussing such systems has been developed in a series of 
papers by Batalin, Fradkin, Fradkina and Tyutin (BFFT) |^, The basic idea is 
to convert, using an iterative prescription, the second class system into a first class 
one by enlarging the original phase space. Once this is achieved it is possible to 
exploit the existing machinery available for quantising first class theories [Qj. The 
use of Dirac brackets is completely avoided. 


The BFFT method, or modified versions of it B B B I, have been used by 
several authors @ to discuss the formulation of specific examples of second class 
theories. These have in turn provided an insight into the general formalism which 
cannot be otherwise gained. But such analyses has been mostly confined to Abelian 
models where the constraints are linear. For nonlinear constraints which occur 
naturally in non-Abelian theories, the situation is rather complicated and explicit 
computations []T| are quite sparse and, at best, limited to a simple group like 
SU{2). The reason for this complication is that, contrary to the example of linear 
constraints, the iterative process in the BFFT approach may not terminate. Nev¬ 
ertheless it was shown by one of us |11] that the nonlinear constraints occurring in 
the nonlinear sigma models could be tackled by a clever choice of the generating 
matrices which forces the iterative process to stop at the first step. Unfortunately, 
such a simplification cannot be done in general |]^]. Then the iterative scheme has 
to be continued indefinitely in which case it is not clear whether consistent solutions 
for all the steps can be obtained from which meaningful closed form expressions for 
the involutive constraints and Hamiltonian may be determined. 


In this paper we systematically develop an algorithm within the BFFT approach 
by which all the iterative corrections necessary to transform the original second class 
nonlinear constraints into strongly involutive constraints can be explicitly computed 
and exhibited in a closed form. We shall discuss this in Section 3 within the context 
of the massive Yang-Mills theory which is a classic example of a second class sys¬ 
tem with nonlinear constraints. The iterative corrections are explicitly computed 
and form an infinite set which is expressed in a closed (exponential-like) form. We 
then use a modified prescription, which is a departure from the conventional BFFT 
ideas, of computing the involutive Hamiltonian. Once again the iterative terms can 
be combined to yield an exponential-like series. The modified prescription men¬ 
tioned above provides great technical simplifications. This has been outlined in 
Section 2 which also contains a short review of the BFFT approach. In our entire 
analysis there is no restriction on either the dimensionality of space time or the 
specific non-Abelian gauge group employed in the model. Since the gauge invari¬ 
ance of the massive Yang-Mills theory is broken by the mass term, it will be clear 
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that our analysis is applicable to all such generic second class systems whose lack of 
gauge invariance is attributed to identical reasons. In Section 4 we first reconsider 
the generalized Stiickelberg formalism [13| of converting the massive Yang-Mills La- 
grangian into a gauge invariant form by the introduction of auxiliary scalars. Next, 
the canonical formalism of the gauge invariant Lagrangian is developed in details. 
Finally, by a series of remarkable algebraic manipulations it is shown that these 
scalars are exactly identified with the additional fields that occur in the correction 
terms found in the BFFT approach. This establishes a one-to-one correspondence 
between the Lagrangian and Hamiltonian embedding prescriptions of converting 
gauge noninvariant into gauge invariant systems. Some appendices are included 
for elaborating calculations relevant for Section 4. Our conclusions are given in 
Section 5. 


2 Brief review of the BFFT formalism 


In this section we present a short review of the Hamiltonian formalism, which will 
be relevant for the subsequent discussion, for converting second class systems into 
true gauge (i.e. first class) systems by enlarging the phase space. Towards the 
end we shall also discuss a deviation from the conventional approach that provides 
considerable technical simplification. 

Let us take a system described by a Hamiltonian Hq in a phase-space with vari¬ 
ables {q^,Pi) where i runs from 1 to N. For simplicity these variables are assumed 
to be bosonic (extension to include fermionic degrees of freedom and to the con¬ 
tinuous case can be done in a straightforward way). It is also supposed that there 
exists second-class constraints only since this is the case that will be investigated. 
Denoting them by Tq, with a = 1,..., M < 2N, we have 

{Ta,Tb} = Aab, (2.1) 


where det(Aab) ^ 0. 

The first objective is to transform these second-class constraints into first-class 
ones. Towards this goal auxiliary variables are introduced, one for each second- 
class constraint (the connection between the number of constraints and the new 
variables in a one-to-one correlation is to keep the same number of the physical 
degrees of freedom in the resulting extended theory), which satisfy a symplectic 
algebra. 


= ( 2 . 2 ) 

where is a constant quantity with det (u;“^) ^ 0. The first class constraints are 
now defined by, 
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Ta = Ta{q,p-,'n ), 


(2.3) 


and satisfy the boundary condition 

Ta{q,p;0) =Ta{q,p) . (2.4) 

A characteristic of these new constraints is that they are assumed to be strongly 
involutive,i-e. 


{Ta, Th} = 0 . 


(2.5) 


The solution of (2.5) can be achieved by considering an expansion of Ta, as 

C50 

fa = Y.T^^\ ( 2 . 6 ) 


n=0 


where is a term of order n in p. Compatibility with the boundary condition (2.4) 
requires that 


= Ta . 


(2.7) 


The replacement of (|2.6| ) into (^) leads to a set of recursive relations, one for each 
coefficient of r/”. We explicitly list the equations for re = 0,1, 2, 


( 2 . 8 ) 

+ {ri 2 ),Tf}(^) + {ri 3 ),rf}(^) = o, (2.10) 


The notations {, }(g,p) and {, }(^) used above represent the parts of the Poisson 
bracket {,} relative to the variables {q,p) and (p). 

The above equations are used iteratively to obtain the corrections (re > 1). 
Equation ( p.sp shall give T^^'>. With this result and (2.9), one calculates T^‘^\ and 
so on. Since T^T ig linear in p we may write 


T^T = Xab{q,p)p'' ■ 


( 2 . 11 ) 


^Sometimes in explicit calculations if no suffix appears in the definition of the Poisson bracket, it will 
imply an evaluation relative to the initial phase space variables {q, p). 
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Introducing this expression into (2.8) and using the boundary condition (^), as 
well as (H and (|2.2| ), we get 


^ab + ^ac Xbd — 0 . 


( 2 . 12 ) 


We notice that this equation contains two unknowns Xab and Usually, first 
of all 0 ;“^ is chosen in such a way that the new variables are unconstrained. It is 
opportune to mention that it is not always possible to make such a choice [14|. In 
consequence, the consistency of the method requires an introduction of other new 
variables in order to transform these constraints also into first-class. This may lead 
to an endless process. However, it is important to emphasize that can be fixed 
anyway. 


After fixing we pass to consider the coefficients Xab- They cannot be ob¬ 
tained unambiguously since, even after fixing a;“^, expression ( ^.12 ) leads to less 
equations than variables. The choice of A’s has therefore to be done in a conve¬ 
nient way [|ll|. 


The knowledge of Xab permits us to obtain Ta^\ If Xab does not depend on 
{q,p), it is easily seen that Ta + Ta^'^ is already strongly involutive and we succeed 
in obtaining Ta- This is what happens for systems with linear constraints. For 
nonlinear constraints, on the other hand, Xab becomes variable dependent which 
necessitates the analysis to be pursued beyond the first iterative step. All the 
subsequent corrections must be explicitly computed, the knowledge of T^\n = 
0,1, 2,...re) leading to the evaluation of from the recursive relations. Once 

again the importance of choosing the proper solution for Xab becomes apparent 
otherwise the series of corrections cannot be put in a closed form and the expression 
for the involutive constraints becomes unintelligible and uninteresting. 


Another point in the Hamiltonian formalism is that any dynamic function A{q,p) 
(for instance, the Hamiltonian) has also to be properly modified in order to be 
strongly involutive with the first-class constraints Ta- Denoting the modified quan¬ 
tity by A{q,p] r/), we then have 


{fa, i} = 0. 


In addition, A has also to satisfy the boundary condition 


^(9,P;0) = A{q,p). 


To obtain A an expansion analogous to (|2.6D is considered, 

OO 

A = J2 > 

n=0 


(2.13) 


(2.14) 


(2.15) 


where A^^') is also a term of order re in r/’s. 
requires that 


Consequently, compatibility with (2.14) 
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= A. 


(2.16) 


The combination of (^), ( 2.13 ) and ( 2.15| ) gives 



(2.17) 



+ {rP,HW}(^) = o, 

(2.18) 





+ {ri3),H«}(^) = o, 

(2.19) 


which correspond to the coefficients of the powers 
expression (2.17|) above gives us 


ry^, ry^, etc., respectively. The 


= -7,<^u;abX^^{T„ A}, 


( 2 . 20 ) 


where ujab and are the inverses of oj°’^ and Xab- 

It was earlier seen that Tq + Ta^'^ was strongly involutive if the coefficients Xab 
do not depend on {q,p). However, the same argument does not necessarily apply in 
this case. Usually we have to calculate other corrections to obtain the final A. Let 
us discuss how this can be systematically done. We consider the general case first. 
The correction A^^'^ comes from equation ( p.l8| ), that we conveniently rewrite as 

= ( 2 . 21 ) 


where 


= {Ta, + {Ti^, + {tP, . (2.22) 


Thus 


= (2.23) 

In the same way, other terms can be obtained. The final general expression reads 

ry“ LOab X^^ G^*") , (2.24) 
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where 


n n—2 

Gi"’ = E {h"-”’, . 

m=0 m=0 

(2.25) 

Although it is possible to convert any dynamical variable in the original phase 
space into its involutive form by the above method, there may be technical problems 
in carrying out this construction particularly if one considers non-Abelian theories. 
In such cases the relevant variable is already quite complicated and this process 
may lead to an arcane structure which would not be illuminating. We suggest the 
following simplification which has been employed, though in a modified form, earlier 
for abelian models |8|. The basic idea is to obtain the involutive forms of the 
initial fields q and p. This can be derived from the previous analysis. Denoting 
these by q and p so that. 


{f,q} = {f,p} = 0. (2.26) 

Now any function of q and p will also be strongly involutive since, 


{T, F{q,p)} 


dP , dP 

0 . 


(2.27) 


Thus if we take any dynamical variable in the original phase space, its involutive 
form can be obtained by the replacement, 

F{q,p) ^ F{q,p) = P{q,p). (2.28) 

It is obvious that the initial boundary condition in the BFFT process, namely, the 
reduction of the involutive function to the original function when the new fields are 
set to zero, remains preserved. 


3 Involutive constraints and Hamiltonian in 
the massive Yang-Mills theory 


The massive Yang-Mills theory is defined by the Lagrangian density, 


L 




(3.1) 
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where the following conventions and notations will be used, 


pa 

flU 

= d,Al-d,Al + gr^^A^^Al, 


[r“, r'*] 

= i , 


tr (r“r'’) 

= 

2 


^rpa'^bc 

= 

(3.2) 


The canonical momentum conjugated to reads, 


TT 


a 


dC 


Hence, ttq is a primary constraint, 


(3.3) 


Tf = vr^ Ri 0 . (3.4) 

In order to look for secondary constraints, we construct the canonical Hamiltonian 
density. 


n, = 


- C, 

^aO 


= ^ 7r“7r“ + 


(3.5) 


The total Hamiltonian [|| is given by, 

Ht = J d^x (7fc + A“7rg) . 


(3.6) 


In Eq. ( |3.6| ) note that the term that appears in (|3.5| ) has been absorbed in 

the A^ttq term by a redefinition of the Lagrange multiplier A“. Consequently, Tic 
occurring in ( |3.6D differs from (^) by this piece. The consistency condition for the 
primary constraint, 


{ttq{x), Ht) 




(3.7) 


yields a secondary constraint, 




= {Dyy + 0 . 


(3.8) 


There are no more constraints since the Poisson algebra ^ of the constraints T“ and 
T 2 is noninvolutive, 


{Tf (x), T^{y)} 
{Tf (x), r|(y)} 
{r“(x), r|(y)} 


0 , 

— 6{x — y ), 

gP^^Dyr^x-y). 


(3.9) 

(3.10) 

(3.11) 


The quantities A^^^(x,y) (^,5 

by, 


1,2) introduced in Eq.(2.1) 


are therefore given 


Ki{x,y) = 0 , 

A?^(x, y) = -m^ - y) = - A^J(x, y ), 

A^^(x, y)=g {Dyf 5{x - y ). (3.12) 


Let us now extend the phase space by introducing the set of new variables 
(t/^“, 7/^“). We consider them as canonical, i.e.. 


{??^“(x), r?^^(y)} = 5“^ 5{x - y ), 

{r/^“(x), y^^{y)} = 0 = {r?^“(x), g^^{y)} . (3.13) 


The symplectic matrix u) defined in Eq. (2.2) therefore has the structure. 



y)- 


(3.14) 


To calculate the first correction of the constraints we have to solve Eq. (2.12), which 
is now written as, 

A^b(x, y) + J dzdz'XTc{x,z)u;<^^<^^{z,z')X%{y,z') = 0. (3.15) 


Using (EH) and ( |3.14| ), three independent relations from ( |3.15| ) are obtained, ^ 


^All algebra will be implemented at equal times. Furthermore, the three dimensional delta function 
S^{x — y) will be written simply as 6{x — y). 

^Detailed reference to the arguments of X and the integral occurring in ( 3.15D will be subsequently 
omitted. 
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A = l,B = l 
A=l, B = 2 
A = 2, B = 2 


0 , ( 3 . 16 ) 

6{x — y), ( 3 - 17 ) 

-gr’’^Dyr5{x-y). ( 3 . 18 ) 

As already discussed the number of equations is less than the number of unknowns. 
Consequently there is an arbitrariness in the solutions even after uj has been fixed. 
At this point it is important to carefully choose the solutions for the variables X so 
that the subsequent algebra is simplified. From the last equation in the above set 
the following choice is practically spelled out, 


X ac \rbc \rbc \rac 

11^12 ~ 

vCLC vbc vbc 

^11^22 “ ^12^21 
VCLC vbc vbc 

^21^22 ~ ^22^21 


= 5^^5{x-z), ( 3 . 19 ) 

^22 = \gf’^{DyY5{x-z). ( 3 . 20 ) 


Using ( 3.16 ), ( 3.17 ) and the above solutions allows us to conclude the results for 
the remaining variables so that the complete X matrix has the structure, 



0 —m? \ , , 

^ab 1 g fabc j 0{X - y) 


( 3 . 21 ) 


Inserting the expressions for the elements of X in Eq. 
the first corrections. 


(|2.11) immediately leads to 


7^(1)“ 

-^2 


gla^^gjabcg2b 


( 3 . 22 ) 


where a compact notation is invoked which will prove highly useful. 


-ac _ j:abCg 2 b 


( 3 . 23 ) 


The next task is to evaluate the second correction. This can be done by using Eq. 


(2.E) which, for the present model, takes the form. 


in. = 0 ■ (3-24) 


Exploiting the expressions for the original constraints (p.4|), ( p.8D and the first cor¬ 
rections ( p.22 ) in the above equation leads to the following conditions for distinct 
values of A and B, 
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A = l,B = l 

+ {Tf= 0 , (3.25) 


A = l,B = 2 

w? {rf\ rf + {Tfr?!' - ^ 5 (Avr*)"}(,) = 0 , (3.26) 


A = 2, B = 2 

mTrr + AS, ^ 

+ {r?'“ + ^ 5 (Avr*)", (A-vr^y + A^} 

+ A'“ + ^5r"(Avr*)",rf'’}(^) 

+ {rf“, A + ^<?^'"(A^')"}(,)=0. 


The above relation (3.27) can be further simplified to yield, 


(3.27) 


{»'“ + is-f"(A tt-)", rf+ {rf 
1 


= -^5V““l?‘'AAvrT5(x-y). 


(3.28) 


It is evident that by choosing the second correction of Tf to be zero achieves consid¬ 
erable simplification. Naturally the same cannot be done for because of (p.28). 


(2)a 


Thus (tl.25D is trivially fulfilled while (|3.26|) will also be satisfied provided Tr 
does not depend on With this hypothesis (p.28 ) is solved leading to the final 
structures, 


T^a ^ Q _ 


(3.29) 

(3.30) 


Note that we have introduced a matrix notation, which will be frequently used, to 
denote the product among the r/ variables. It is now useful to explicitly write the 
involutive constraints obtained up to the second iterative step. These are given by, 


A = 7rg-mA^“ + --- (3.31) 

A = m2AS + (AvrT+r?'“ + ^5r'=(A7r*)^ 

+ ^(fr(AvrT + --- (3.32) 
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It is already seen that while T“ has a simple structure, the correction terms in 
T 2 generate an exponential-like series. This is further elaborated by considering the 
third step of the method. The equation we have to solve occurs in ( p.iop which, in 
the present instance, reads, 

{T%, rf+ {Tf T^} + 

+ {rf “, + {Tf\ = 0 . (3.33) 

We thus have the following relations for distinct values of A and B, 


A = 1 R = 1 




(3.34) 


A = l, B = 2 

^2 r 2a t ,^2 rT-i(3)a 16 


m-W\ Tr} + m^{Tr\ + - 5 ^ (AvrT}(,) = 0 , (3.35) 


A = 2, B = 2 


+ 2 5 r" (Avr*)", rf + {rr“, 7?- + - g r (A7r‘)'^}(,) 


■’(3)^^ ^16 I ^ /■ D 


^3 jbcd _j_ _ JO'Cd 

+ ^ {Dyr 6{x - y ). (3.36) 


Using our algorithm of taking a vanishing correction for Tf, it is simple to observe 
that these equations will be satisfied by choosing. 


T(3)a ^ Q ^ 

3 

rf (fj3ff (BiTrY . (3.37) 


This iterative process can be extended to arbitrary orders and the hnal expressions 
for the involutive constraints are given by, 


-‘I ~ ”o m 7] , 

oo 

Al + y^^ + Y^ 


n=0 


iavY 


-I ab 


(n + 1)!_ 


{Di7r‘ 


i\b 


(3.38) 
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This completes the first part of the analysis. It is now necessary to construct 
the involutive Hamiltonian. We adopt our modified prescription of first doing this 
construction for those phase space variables occurring in the initial canonical Hamil¬ 
tonian. The calculation for vrf is now given here in details. The first step is to com¬ 
pute the inverses of the matrices ( |3.14 ) and (3.21), denoted by and 

respectively. These are given by, 


oM - 9) 


( 3 . 39 ) 




— I 2 m,' 




1 mb 


- -L(5' 


0 


5{x - y) 


( 3 . 40 ) 


Using the elements of these matrices the first correction is obtained from (2.20), 


^ /^i 0 )ba I 9 fbcd ^ 2 b t „^^c/'^(O)rfa ^ 2 b/^i 0 )ba 


( 3 . 41 ) 


where the generating functions are defined by. 


G 

G 


Sf“ = {r^7r“}=0, 

ma ^ {r|,vr“}=5/'“‘'vrf, 


( 3 . 42 ) 

( 3 . 43 ) 


Inserting these in Eq. (3.41) the explicit result for the first correction follows. 




( 3 . 44 ) 


This will be used to derive the next correction. The result for an arbitrary iteration 


is given from the general formula ( 2.24 ), which is further simplified by inserting the 
matrices ( 3.3£ ) and ( 3.40|) 


vr. 


{n+l)a _ 


^ „Ab , .be vBcCd r^{ri)da 

^ Mn)ba _ 1 

/ -I N O'/ 1 I 

(n -|- l)m^ 


—be f j-i i\b ^(ji)ea 


2{n + 1) m? 




^ „2f) Mn)ba 

'/ '^2 


n + l 


( 3 . 45 ) 


The structure of the second correction is therefore given by. 


( 2 )a _ 1 16 ^(l) 6 a 9 6 c / m i ',6 ^(l)ca 1 26 

(A7r)Gi, --7? Gs, , 


( 3 . 46 ) 
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where the generating functions, defined in Eq. (p.25), are, 


^{l)ha 

^{l)ha 

'^2i 


0, 

(3.47) 

{r|, vrf)“} + <} + {rfvrf )“}(^) , 


g2 jbcdjaec ^2e^d ^ 9 _ jbcd jdah ^2c^h _ 

(3.48) 


Using these results in Eq. (|3.4(:), we find the explicit form of the second correction. 


vr. 




(3.49) 


Similarly, the expression for the third correction looks like. 



1 

2>m? 


lb /~i{‘i)ba 


+ 


9 

Qm? 




a 


{2)ca 


li 


^ 2b 

3^ ^2* , 


(3.50) 


where. 


^(2)ba 

^[2)ba 

^2% 


<} + 7rW“} + {^2^ vrf 

„3 

jbcd^aef^dgh ^2c^2e 


jacdjdefjbgh ^2c^2e 

2 


(3.51) 


(3.52) 


Using the above results in (3.5C|), we find. 




(3.53) 


It is clear that the structure for the corrections to the constraints has considerably 
simplified the algebra. There is no correction to the first generating function while 
the second acquires a typical form. This leads to a correction in vr? which is easily 
generalized to arbitrary orders. The final result for the involutive expression is 
thereby given by. 
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(3.54) 


OO 



Likewise it is possible to compute the involutive expressions for the other variables. 
Since same steps are employed which lead to identical simplifications, the details 
are omitted. The final expression for the involutive is given by, 


n=l 

OO 1 

.a(n + l)! 

The corresponding compntation for is slightly tricky since the normalisation is 
no longer given by a simple factorial. We quote the final involutive structure. 


OO 

la _ 4a I ^ 4 Wa 
n=l 

1 OO 

= + — E - 9" bn {Dyr 

n=l 

where the coefficients are found from the recursive relation, 


(3.56) 


hi 


1 

n 


'n-l + 


1 n 
(n + 1)!- 


1 

2 ■ 


n > 1, 


(3.57) 


The involntive Hamiltonian H is now obtained directly from the canonical Hamil¬ 
tonian that appears in Eq. (|3.6|) by snbstituting the initial phase space variables by 
their corresponding involutive expressions, i.e.. 


H{AQ,Ai,^i) 


Hq^Aq^ Ai^ TTi) 


^0 ^ ^0 

Ai > Ai 

TTi > TT^ 


(3.58) 


This completes the BFFT convertion of the second class massive Yang-Mills theory 
into a true (first-class) gauge theory. The constraints (|3.38| ) and Hamiltonian ( ^.58 ) 
satisfy a strongly involutive algebra in the extended phase space. 
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4 Generalized Stiickelberg formulation 


It is well known that a gauge noninvariant Lagrangian can be converted into a 
gauge invariant form by introducing auxiliary scalars. This is the Stiickelberg |]T] 


formalism by which it is possible to discuss a gauge invariant formulation of the 
massive Maxwell theory, usually called the Proca model. A non-Abelian generalisa¬ 
tion of the Stiickelberg mechanism was first proposed by Kunimasa and Goto 
It was subsequently used to analyse a gauge invariant Lagrangian formulation of 
the massive Yang-Mills theory by a number of authors 0,0ii- These analyses 
are presented either from a geometric or path integral viewpoint. In this section 
we shall first develop in details, using the generalized Stiickelberg prescription, the 
canonical formalism of the gauge invariant formulation of the massive Yang-Mills 
theory in the coordinate language. Subsequently it will be shown that the auxiliary 
scalars (and their conjugates) introduced in this case are just the canonical pairs 
(^ia,^ 2 a) yggjj BFFT approach. 

The pure (massless) Yang-Mills Lagrangian is invariant under the gauge trans¬ 
formations, 


A? 




where 9°'{x) are the group parameters and, 

^“( 0 ) = V^\9)d^6K 

The transformation matrices U and V are given by. 


C/“^(0) 


V^^(9) 


- 1 ab 

exp(- g9) 


E 

n=0 


(n -h I)! 


ab 


where. 


^ab _ jaco q 


acb nc 


(4.1) 


(4.2) 


(4.3) 

(4.4) 


(4.5) 


uses a notation exploited previously in ( 3.23 ). From the above relations it is simple 
to deduce the result for infinitesimal variations. 


6Af, = U{9)6A^ + D^ [A®] {V{9) S 9) 


(4.6) 


and Dfj_ is the covariant derivative already introduced in (^(^). The inclusion of a 
mass term in the pure Yang-Mills Lagrangian, as has been done in Eq. (^), breaks 
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the above gauge symmetry. Nevertheless, it is straightforward to incorporate a 
non-Abelian gauge invariance into the massive Yang-Mills Lagrangian by extending 
the configuration space. This is the content of the generalized Stiickelberg mecha¬ 
nism. Using (4.2) and ( [4.6| ) it can be shown [17| that by construction, the following 
Lagrangian, 


1 ?77^ 

L = -- ^ ( 4 . 7 ) 

is invariant under the infinitesimal non-Abelian gauge transformations, 

Qa —^ r + A“, (4.8) 


where A is the gauge parameter. 


Let us next develop the canonical formalism for the gauge invariant Lagrangian 
(4.7). For subsequent calculations it is useful to introduce the matrix W°'^{6) where. 


w^\e) 


- [ {-gOr 


(4.9) 


Up to a four divergence the Lagrangian (4.7) may be written as, 


L = 

- ^ Fl, + ^(71^+ B]) - vf? A, [df^9 + W{9) df^9] , 


= 

- 1 Flu + ^ (tIm + fD - m2 A, W(9)8^9 + m2 9 d,A^ . 

(4.10) 

The canonical 

momenta are. 



dC 



< = =m^9^, 

0 9^0,a 

(4.11) 


dC 



TPCi 

(4.12) 


0 

TT^ = ^ = m2 9'^ - m2 A[1. 

(4.13) 


There is one primary constraint, 

Tf = TT^ - m2 (9“ 0 . (4.14) 
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The above manipulations of isolating the identity component from V°‘^{9) by 
defining W°'^{9) (|4.9| ) and then expressing the Lagrangian up to a four divergence 
now become clear. Our motivation was to obtain a constraint that could be iden- 


tihed with Tf in (3.38). Unless this manipulation was done the primary constraint 


obtained from (4.7) would just be ttq ~ 0 and the mapping fails. 

For the computation of the secondary constraint, it is necessary to obtain the 
expression for the canonical Hamiltonian. This is given by, 


H = 


J (fx ^TTq + TTj + TTg 0“ — cj 


m 


= Hc + ^ J (fx 9‘^9‘^ + AH , 


(4.15) 


where He is the canonical Hamiltonian that appears in Eq. (| 
massive Yang-Mills theory, and, 


AH = mr 


'Sx[]-B^ 




for the usual 


(4.16) 


involves terms not depending on either Aq or 9 so that it does not influence the 
computation of the secondary constraint. In order to simplify the velocity dependent 


term in ( 4.15 ) the first step is to invert ( 4.13 ) so that the velocity 0“ is expressed 
in terms of the momenta tt^. This is done in Appendix A. Next, using the results 
of this Appendix, the desired simplification is done in Appendix B. Using (B^) in 


(4.15) the final form of the Hamiltonian is obtained. 


H = Hc + AH + ^ I d^xT^T^^{-lY (A") 


ga 


(4.17) 


n=0 


Time conserving the primary constraint yields the secondary constraint T 2 , 

{Tf, ft} = - rrA 0“, Ft} = 0 , (4.18) 


where. 


+ Ar“ 


(4.19) 


involves the sum of the secondary constraint (3.8) occurring in the usual massive 
Yang-Mills theory and an extra piece. 


A ^da , 'Z Ad (-\(la 

i\i = — 'Kq — TTq X + ^ ^0 ^ ^ 


^2 Ad r\da 


(4.20) 


The structures x and have been defined as well as simplified in Appendix C (C.15 
and p.26|) . 
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We now make certain observations. An inspection shows that the constraint 


(4.14) gets mapped on to the first constraint in (3.38) by a simple identification of 
fields. However the same thing cannot be done for the other constraint. Neverthe¬ 
less, as we shall now show by a series of manipulations, a correspondence between 
these constraints can also be achieved. Using (p.26|) in (4.20|) we obtain. 


Ar“ = 0'^ - TT^ , 


(4.21) 


where 


= + (4.22) 

0n the constraint surface T 2 = 0 we find the solution,^ 

d ^ 

0“ = - (Avt*) ^(-l)'^ , 

n=0 

= - + (4.23) 

Defining the inverse by I so that, 

^ jbc ^ ^ac _ ^4 24 ) 

Expressing this as a series, 

00 

= + [(50)1“' (4.25) 

n=l 


and inserting it in (4.24) leads to the following condition for f3n] 


E [(50)1“' + ^2“' + E [(5011“' = 0 


n=l 


n=l 


(4.26) 


From the above relations, therefore, the extra piece AT“ simplifies to 


AA = - (Avr*) “ - Tie , 


LXJ 

(Avr*)“E/^- [(50l1““-<- 


n=0 


(4.27) 


"‘it is crucial to note that the equality between y and U is essential to derive this result, which plays 
a central role in the forthcoming analysis. 
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The final task is to compute fdn by solving Eq. (4.21;). This is done in Appendix D. 


Using Eq. (D.4) the final expression for the constraints is given, 




T id 
2 


la ,^2/ia 


= - m 

= vtq — , 


= n- 


+ (AvT*) , 


n=l 


m^A^-TTg + Y^ an [{gO)Y'^ (Att*) , 


n=0 

oo 


= m 




n=0 - 




“I a6 


(n + 1)! 


(Avt*)' , 


(4.28) 


where (A.2) has been used to obtain the last line. If we make the following identi¬ 
fications with the fields introduced in the BEET approach, 


e° 


n 


2a 


Til ^ 


■ 7 ] 


la 


(4.29) 


the canonical algebra is preserved, 

{^“(x), 4(y)} = - {??^“(x), r]^\y)} = <5“'’ 6{x - y). (4.30) 


Now the above results for Tf and are exactly identifiable with the corresponding 
expressions given in ( |3.38| ). This shows that the auxiliary scalars in the generalized 
Stfickelberg formalism are exactly mapped on to the BEET fields. 


Before closing this section it may be worthwhile to point out that the constraint 
the form ( 4.2^ ) (and not as it occurs in ( 4.2C| )) reveals that it is the generator of 
gauge transformations because. 


m 


j d^x{f^{x)\-{x),e\y)} = X\y), 

Jd^x{fi{x)X‘^{x),A'{{y)} = (^"“A“(y)) , (4.31) 


which correctly reproduces (|4.8|). 
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5 Conclusion 


We have discussed a systematic method, within the Batalin, Pradkin, Pradkina and 
Tyutin (BPPT) [^, ^ approach, of converting the massive Yang- Mills theory to 
a gauge invariant theory by extending the phase space. Exploiting an intelligent 
choice for the symplectic matrix and the generating matrix the infinite 
number of iterative corrections necessary for obtaining the strongly involutive con¬ 
straints were considerably simplified. These corrections were explicitly computed 
and expressed in a closed (exponential-like) form. In obtaining the strongly in¬ 
volutive Hamiltonian, on the other hand, the conventional BPPT approach was 
modified. Pirst, the strongly involutive forms for the initial phase space variables 
were deduced. Once again infinite sets of iterative corrections were necessary which 
were computed and put in exponential-like series. Then the canonical Hamiltonian 
of the massive Yang-Mills theory was rewritten, replacing the original phase space 
variables by their corresponding involutive expressions. This directly gave the cher¬ 
ished form of the involutive Hamiltonian. Indeed any dynamical variable in the 
original theory can be converted into its involutive form by this technique. 


In the latter half of the paper the generalized Stiickelberg formalism |13] of 
changing the gauge noninvariant massive Yang-Mills Lagrangian into a gauge invari¬ 
ant form by introducing auxiliary scalars was reconsidered. The complete canonical 
formalism of the gauge invariant Lagrangian was developed in the coordinate basis. 
The explicit structures of the involutive constraints and Hamiltonian were deter¬ 
mined. Subsequently it was shown by a series of algebraic simplifications that the 
auxiliary Stiickelberg scalars and their conjugates were exactly identified with the 
additional canonical pairs of fields, defined in the extended phase space, invoked ear¬ 
lier in the context of the BPPT analysis. By this identification a mapping between 
the Lagrangian and Hamiltonian embedding prescriptions of transforming gauge 
noninvariant into gauge invariant systems based on the generalized Stiickelberg and 
the BPPT approaches, respectively, was established. It should be stressed that this 
mapping is independent of either the space-time dimensionality or the specific non- 
Abelian gauge group employed in the analysis. In this sense, therefore, the present 
paper extended and generalized similar correspondences reported earlier i 0] for 
abelian groups. In this instance the abelian result follows trivially by setting the 
structure constants to zero. 


It is clear that just as the generalisation of the usual Stiickelberg formalism 
flstl from Abelian to non-Abelian theories is nontrivial, the same is true in the 
BPPT formalism. Indeed a distinctive feature of transforming non-Abelian second 
class into first class systems, in contrast to abelian theories, was that all orders of 
iteration in the BPPT approach were mandatory. But our approach provided an 
algorithm of systematically computing these corrections, and expressing them in 
closed forms. An interesting application of the method developed here would be to 
study the bosonisation and duality among non-Abelian theories in higher dimensions 
where the conventional master Lagrangian approach |19] had failed whereas, that 
based on the generalised Stiickelberg formalism was effective |M|. 
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As final remarks we mention that the Hamiltonian formalism for the massive 
Yang-Mills theory, regarded as a second class system, was originally developed 
in [^. The corresponding first class (canonical) interpretation in terms of the 
Stfickelberg approach was earlier presented in [l^ which, contrary to the present 
work, starts from an embedded Lagrangian and does not use the elaborate BFFT 
method. Incidentally, the idea of converting second class systems their first class 
forms by extending the phase space, which is the crux of this method, was initially 
suggested in |]^ . 
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Appendix A 


This Appendix is devoted to invert ( 4.13| ) so that the velocities are expressed in 
terms of the momenta. The product of the V matrices ^ is simplified by using (|4.4D , 




QO 


E 

n=0 


1 

(n -|- l)(2n -|- 1)! ’ 


V^^{9) 


P^“(0), 


(A.l) 


where the symmetry follows from the definition of 6 given in Eq. (4.5) 


since, 


in 


ab 


(_l)n 


(A.2) 


It is once again useful to separate the n = 0 contribution. 




A“'’(0) 


“ r (ffg)^» 

[(n + l)(2n + 1)! 


A^“(0). 


(A.3) 

(A.4) 


Defining 

®In a geometric language this product has the meaning of a two sided invariant metric on the group 
of rotations [l^ . 
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(A.5) 




and using the above relations enables Eq. (|4.13 ) to be expressed as, 

r“ = r+ 


The solution for 0“ is now given by, 

OO 

0“ = r^^(-l)” (An)ba ^ 
n=0 


which may be verified by directly substituting in the RHS of (A.t), 


OO OO ^ 

pfe ^ (_ i)n ^ ^ “ _ 

n=0 n=0 


Changing the sum in the second term from n to m = n + 1 leads to 

OO OO 

^prn\^ha _ pa 

n=0 m=l 


(A.6) 


(A.7) 


(A.8) 


(A.9) 


which reproduces the LHS of (A.6). 


Appendix B 


The velocity dependent factor in the Hamiltonian (4.15) will be expressed here in 
terms of the phase space variables. This factor is written as 


ycuycbgaQb ^ (^ab ^ ^ab\ gagb ^ 


OO OO 

= ( 5 ^^ + r ^ (-1)” (A”)“ T'^ ^ (-1)”^ (A”^)'^'’, (B.l) 


n=0 


m=0 


where we have used (A.l), (A.3) and (A.7). The RHS simplifies to. 


OO OO 

r ^(-l)*" (A^)^^ T'^ ^ (-I)’” (A”^)'^^ 

n=0 m=0 

00 ^00 

+r^^(-i)'^ (a"+^)'' t'^ ^(-i)'" (A'")'^'’. 

n=0 ?n=0 


(B.2) 
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Consider the second term. Replace the first sum from n to = n + 1, 

oo ^ oo 

- r ^ (-1)™' (A™y ^ (-i)’” . (B.3) 

m'=l m=0 

Adding this with the first term in Eq. (|B.2| ) yields the final result 

oo 

ycaycbQaQb ^ papfe ^ (-1)™ (A™)“'' . (B.4) 

m=0 


Appendix C 


The matrices x ^md fl appearing in the constraint 
details. We first treat x whereas Q is taken up from 
elements of x are given by, 

OO OO 

^ab ^ ^ y,ac ^ ^ 

n=l n=0 

where W and A have been introduced in ( [4.S| ) and ( |A.4D , respectively. Consider the 
expression, 


(4.19) are considered here in 
(C.17) onwards. The matrix 


OO 

^ (-1)" (A")“'' = (5“'^ + A'^'') “ . (C.2) 

n=0 


Writing this as a series in 9, 

OO OO 

£ (_1)- = £ a [(<70)”]“' , (C.3) 

n=0 n=0 


it is possible to determine Cn using (|C.2[) , 


n=0 


(C.4) 


It is easy to verify that only even powers of {gO) contribute. The final result is 


^ (-1)^ (A’^)“'’ 


n=0 


oo 

n=0 


\sS) 


(C.5) 


where 
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Co = 1, 


Cn — ^ ^ Cn—m ^ 0 ) , 


m=l 


Am. — 


{m + l)(2m + 1)! 


(C.6) 


We now simplify the second term in (C.l); 


n=0 p=l n=0 


2n 


c6 


(C.7) 


with. 


(-1)" 

(P+ 1)! 


(C.8) 


Let ns express the product occurring in the RHS of (C^) as a power series where 
the odd and even terms are explicitly separated, 


gO 


ab oo 


2 ! 


n=l 


ab 


+ J2dn (ge)^^ +Y.fn W 


n=l 


\2n-\-l 


ab 


(C.9) 


A straightforward comparison of terms fixes dn and fn to be, 


n—1 


dn — ^ ' 0.2{n— 


m) Cm ) 


m=0 


fn = E ®2(n—m)+l Cm 

m=0 


(C.IO) 


Writing the complete expansion for /„ yields. 


fn — ®2n+l Cq “t“ 0.271—1 Cl + • ■ ■ + Hi Cn , 

= 02n+l Co + 0271-1 Cl + • • • — ai (Ai Cn-1 + ' • • + An Cq) , (C.ll) 


where we have inserted the expansion of Cn in the last term. Now using the identity 

Ol An — 0.271+1 7 
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that follows from the respective definitions (C.6) and (C.8), it is seen that the last 
factor in the parenthesis cancels the first term in /„. Likewise, there will be a pair 
wise cancellations among all terms so that, 


fn = 0. 


Collecting all pieces together, we obtain, 

OO n 

x"” = (A”)“* - ^ + E <4. [(9«) 


n=l 


ba OO 

E 

n=l 


2n 


ba 


(C.13) 


(C.14) 


Using the result (C.5) and the symmetry properties of is further simplified. 


.,ab _ 9 nab 


where. 


+ E (5^) 


n=l 


\2n 


ab 


(C.15) 


UJn — Cfi + dfi , 
n 

^ ^ ®2(n—m) Cm , 
m=0 


(C.16) 


which follows from the explicit expression for dn given in 


(IC.IO ) and oq = 1. 


We now consider the matrix which is defined by. 


n=0 n=0 


_ _ ^^o.c J gbc _ I \ ^ , . r/'^/3\2ni^c 


(C.17) 


n=l 


where the second line follows from the previous analysis in this Appendix. Putting 
the explicit form for W leads to, 


00 1 00 

Si"* = -{E“p[(s«)’’l“‘ + 5E“pl(9«Wl“* 

P=1 p=l 

OO OO 

+ E “P E “■> [(9«^"”l‘*} ■ (C-18) 

p=l n=l 


The hrst two terms are combined by a simple change in the summation p —s- p+1. 
It leads to, 
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oo oo 

p=l n=l 



q 

2iq + 2) 


a, [{gey+Y + 



(C.19) 


We now prove that terms involving odd powers of g in the two summations con¬ 
spire to vanish. An arbitrary odd powered term {g6y"'~^^ (n = 1 , 2 , ...) has the 
coefficient (apart from a minus sign), 


P2n+1 — ®2n-l + 02n-3 ^2 + ' ‘ ‘ + ^3 ^n-1 + i^n + 


n 

2n + 2 


a2n ■ 


(C.20) 


Inserting the expansion for ujn, given in (p.l 6 |), this yields, 


P2n+l 


a-2n-l (a2Co + CLqCi) + a2n-3 (a4Co + a2Ci + aoC2) 

Ti 

+ •••-1-01 (a2nCo + fl2n-2Cl + • • • + O-oCn) + 7 ;- CL2n ■ (C.21) 

2 n + 2 


The term oiOoCn in the last parenthesis is further elaborated by decomposing Cn, 
OlOQCn — (AiCji_i + A.2Cn—2 + • • • + AjiCq) . (C.22) 


Using the identity ( p.l2| ), it is found that the penultimate term in (p.22|) cancels 
the second factor in the first parenthesis of (p.21j) . Indeed, just as discussed in 
Appendix B, there will be pair wise cancellation except for the last factor in (C.22) 
and the first one in the final parenthesis of ( C.21 ). Combining these remaining 
pieces. 


P2n+1 — a2n Cq — dO Cq + 77 -—- a2n — 0 , (C.23) 

2 n + 2 

which follows from the explicit expressions. Hence, as announced before, terms 
with odd powers of g vanish. Next, looking at even powered terms, the coefficient 
of O (< 7 ^"') (n = 1 , 2 , ...) is given by, 


P2n — — {“2n-2 d>i + a2n-4 ^2 + ' • • + 04 O'n-2 + O2 ^n-l} 


2 n- 1 

2 ( 2 n + l) ■ 

(C.24) 


By following analogous steps it may be shown R that this simplifies to 


®We shall subsequently discuss in Appendix D ( D.8 to D.ll ) how this result also follows from a different 
viewpoint. 
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(C.25) 


n 

P2n — 'y ' 0,2(n—m) ■ 

m=0 


Hence we find the remarkable result, 





n=l 


X 


ab 


(C.26) 


Appendix D 


The explicit solution for j3n will be obtained by first expressing (4.26|) in full form 


X! [{gGT] + + '^^n {gO) 


n=l 


2! 1 
n=l 

gQac 


\2n 


ab 


+ ^+E ‘-n \{g9rr E /5n mr' = o. (d.i) 


n=l 


n=l 


By equating coefficients of the first few terms, we find, 

/^l + ^ = 0 5 

/32 + wi + ^ = 0 , 

02 

0Z + + ^101 = 0, 

0-\ 

0A + UJ2 0- 0- ^102 = 0 , 


(D.2) 


The first equation gives 0i which is used to obtain 02 from the next equation. 
This iterative process can be continued to obtain all the /3’s. Indeed an explicit 
computation, using the value of LOn from (C.16), shows 
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/3i 

P 2 

/?3 

(3a 


1 

~2\' 

1 

1 

“1’ 

1 


(D.3) 


which suggests a general solution for j3n\ 

(-ir 


Pn = 


(n + 1)! 


— Q-r; 


(D.4) 


It is now simple to prove this result explicitly. Take the coefficient of the 
term (n = 1, 2, • • •), 


P2n+1 + ^ + 


^ ^ P2{n—m)+l — ^ . 

m=l 


(D.5) 


Setting the solution (p.4D for /3n, the LHS of the above equation yields, 


, “271 , 

“2n+l + a2(n-m)+l 


m=l 


— 1 0,2n+l + 


02n 


+ 0-2n-l + ^2 fl2n-3 + ' ‘ ‘ + <^n ■ 


(D.6) 


The series separated from the parenthesis has already been evaluated in Appendix 
C (see p.2(]| and C.23). Inserting this result in (p.6|) we obtain, 


^2n ^ ^ 

(X2tiA-\ \ - 0^271 — U , 

+ 2 2n + 2 


(D.7) 


obtained by using (^). Thus the consistency of the solution for Pn is verified. 

It is instructive to observe the coefficients of the 0{g^^) term (n = 2, 3, ...) 
which is given by, 


P2n + U}n+ ^ UJm P2{n-m) — 0 • 


n—1 


(D.8) 


m=l 


Inserting the solution for /?„, it follows. 
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(D.9) 


n—1 

^ ^2{n—m) ®2n 

m=l 


Q2n-1 

2 ! 


— Wn • 


It may be easily checked that, 


02n + 


02n-l 

2 ! 


2n- 1 
2(2n + l) 


(D.IO) 


and we obtain, 


n—1 

~ ^ ^ 0‘2{n—m) 

m=l 


2n- 1 
2(2n + l) 


(D.ll) 


which just reproduces the result ( C.24 ) and ( p.25|) Q. 
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